home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlarfg.z / dlarfg
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAARRRRFFFFGGGG((((3333FFFF))))                                                          DDDDLLLLAAAARRRRFFFFGGGG((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLARFG - generate a real elementary reflector H of order n, such that   H
  10.      * ( alpha ) = ( beta ), H' * H = I
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLARFG( N, ALPHA, X, INCX, TAU )
  14.  
  15.          INTEGER        INCX, N
  16.  
  17.          DOUBLE         PRECISION ALPHA, TAU
  18.  
  19.          DOUBLE         PRECISION X( * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      DLARFG generates a real elementary reflector H of order n, such that
  23.                (   x   )   (   0  )
  24.  
  25.      where alpha and beta are scalars, and x is an (n-1)-element real vector.
  26.      H is represented in the form
  27.  
  28.            H = I - tau * ( 1 ) * ( 1 v' ) ,
  29.                          ( v )
  30.  
  31.      where tau is a real scalar and v is a real (n-1)-element
  32.      vector.
  33.  
  34.      If the elements of x are all zero, then tau = 0 and H is taken to be the
  35.      unit matrix.
  36.  
  37.      Otherwise  1 <= tau <= 2.
  38.  
  39.  
  40. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  41.      N       (input) INTEGER
  42.              The order of the elementary reflector.
  43.  
  44.      ALPHA   (input/output) DOUBLE PRECISION
  45.              On entry, the value alpha.  On exit, it is overwritten with the
  46.              value beta.
  47.  
  48.      X       (input/output) DOUBLE PRECISION array, dimension
  49.              (1+(N-2)*abs(INCX)) On entry, the vector x.  On exit, it is
  50.              overwritten with the vector v.
  51.  
  52.      INCX    (input) INTEGER
  53.              The increment between elements of X. INCX > 0.
  54.  
  55.      TAU     (output) DOUBLE PRECISION
  56.              The value tau.
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.